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ABSTRACT 

We present an abundance analysis of six main sequence turnofF, subgiant, and giant branch stars toward the 
Galactic bulge that were observed with Keck/HIRES during microlensing events. This is an early look at 
the first detailed chemical analysis of main sequence stars in the Galactic bulge. Lensing events allow the 
effective aperture of Keck to be increased beyond its current dimensions; although, some events still stretched 
its spectroscopic capabilities. Future large telescopes with high resolution and high throughput spectrometers 
will allow the study of abundances in distant stellar populations and in less evolved stars with greater ease. 
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1. INTRODUCTION 

The distance to the Galactic bulge has until recently kept the prospects of observing directly the chemical 
properties of its constituent stars just past the edge of possibility. Even now, with large-aperture telescopes 
such as the Kecks, the VLT, the HET, and the Gemini twins among others, only the brightest giants in the Bulge 
have been observable.^ A novel approach employed by the MACHO microlensing team has taken advantage of 
the temporary brightening of the source star during the lensing event to effectively decrease that once prohibitive 
distance to the Galactic center and allow for the observation of less evolved stars. For example, Minniti et al.^ 
first published results in 1997 on the presence of a Li I feature observed in a hot {Tcs = 6000 K) main sequence 
turn-off star. In this paper we discuss the full analysis of the that spectrum along with the spectra of five 
other stars observed during lensing events between 1997 and 1999. The results to date are still preliminary; 
however, they offer a glimpse into the diverse chemical evolution scenarios that can be found in and around the 
center of the Milky Way. Future telescopes with effective apertures in the 30-m or more range will make such 
observations commonplace rather than rare. 

2. OBSERVATIONS 

Between August 1997 and July 1999 we observed with the Keck I telescope and the HIRES echelle spectrometer 
18 microlensing events that occurred when a small compact object (e.g., a normal star) crossed the line of sight 
between earth and a Bulge member, causing a temporary increase in the brightness of the latter by gravitational 
lensing. Of these, only six spectra had a sufficiently high signal-to-noise ratio to undergo a complete analysis 
that allowed us to determine: 1) radial velocities; 2) atmospheric parameters; and 3) chemical abundances. 
These six stars are described in Table 1, where they are identified by their lensing event number and their 
MAGHO id in the first two columns, their unlensed magnitude and color in the next two columns, the epoch of 
their observation and the exposure time in the fifth and sixth columns, and an estimate of the signal-to-noise 
ratio in the last column. The (V — R)o colors are corrected based on recent extinction maps by Popowski, 
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Cook, and Becker.^ We note that the events 97-BLG-45 and 97-BLG-47 were indeed the same star identified 
separately in overlapping MACHO fields. The event 98-BLG-6 was observed spectroscopically during two 
separate epochs in 1998, with the spectra being combined after ensuring no variation in the relative quality 
of each set of observations. The signal-to-noise ratio of this event represents that of the combined (16,200 s) 
spectra. Information regarding these and other microlensing events can be found on the MACHO project home 
page at http://www.macho.mcinaster.ca/. 



Table 1. Program star data 



Lensing Event 


F.T.S 


V 


{V - R)o 


Epoch 


t'cxp (^) 


S/N 


97-BLG-45/47 


142.27650.6057/136.27650.2370 


19.8 


0.42 


28 Aug 1997 


10800 


40 


97-BLG-56 


403.47671.57 


16.1 


0.85 


27 Aug 1997 


450 


80 


98-BLG-6 


402.48103.1719 


19.5 


0.59 


20 June 1998 


10800 


60 










18 Aug 1998 


5400 




99-BLG-l 


121.22423.1032 


18.8 


0.58 


6 July 1999 


7200 


65 


99-BLG-8 


403.47849.756 


15.8 


0.82 


6 July 1999 


1800 


130 


99-BLG-22 


109.20893.3423 


19.6 


0.79 


7 July 1999 


2700 


90 



Our setup provided us with wavelength coverage from 4780 A to 7150 A with gaps between the orders 
that increased with wavelength. This allowed for the observation of many atomic absorption lines; however, it 
missed the important [O I] feature at A6300. Wc binned the data 2 x 2 in order to increase signal and decrease 
readout times. Combined with typical seeing of 1", this produced a resolution of R ~ 29,000; sufiicient for doing 
abundances work. 

The data were reduced using the batch mode of MAKEE, written by T. Barlow and modified by him to 
accommodate our 2x2 binning format. MAKEE takes raw images as input and, after correcting for bias and 
flat fielding, extracts and wavelength calibrates the spectra. We then converted the MAKEE fits output files 
into an IRAF readable format using the MAKEE routine linear. Once the spectra were extracted we used the 
spectrum of a rapidly rotating hot star to correct for telluric lines by dividing the object spectrum by that of 
the rapid rotator, with the two spectra scaled and shifted by hand in the IRAF routine telluric. This removed 
most of the blaze function, allowing us to fit a low order polynomial to the data to produce flat, normalized 
spectra. Wc then measured the equivalent widths of lines by fitting gaussians with the IRAF routine splot. To 
test the quality of the data reduction, we reduced the spectra of many comparison stars that were taken with 
our program stars and that have equivalent widths published in the literature (e.g., Edvardsson et al.*). We 
found the agreement between our equivalent widths and the published data to be generally excellent, as shown 
in Fig. 1 for a sample star. 

3. ANALYSIS 

We began our analysis by determining heliocentric-corrected radial velocities. To do this we used the IRAF 
routine rvidlines and the iron line list of Nave et al.,^ from which we selected approximately 450 Fe I lines 
with accurate wavelengths that fell within our spectral coverage. During each epoch of observations we observed 
comparison stars with well determined radial velocities, and our analysis of the data showed excellent agreement 
between the radial velocities determined for these and their published values. The final column of Table 2 lists 
the radial velocities of our program stars. 

Using the stellar atmosphere code MARCS (Gustafsson et al.^) and the LTE abundance analysis code MOOG 
(Sneden^), we determined model atmosphere parameters for our program stars by using the Fe I and Fe II lines 
to achieve: 1) no trend in abundance versus excitation potential to determine T^g; 2) no trend in abundance 
versus reduced width to determine microturbulence (^); and 3) agreement in abundances from neutral and 
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Figure 1. Equivalent width comparison between our data and Edvardsson et al.* for the star HR 5019. The dotted line 
represents perfect agreement. 

ionized atoms to determine log g. Wc used ^ 49 Fe I lines from Nave ct al.^ and ^ 14 Fe II lines from Bicmont 
et al.,^ with the latter's \oggf values corrected to loge(Fe) = 7.52 (a scale relative to hydrogen where log(H) is 
taken to be 12.00). The estimated uncertainties are ±200 K in Tgff, ±0.5 dex in log g, and ±0.2 km in ^. 
The derived atmospheric parameters are listed in Table 2. 

We could have used the V — R colors and the scales of Alonso et al.^ or Houdashelt, Bell, and Sweigart^° 
plus some assumptions about the distance modulus of each star to determine atmospheric parameters, but the 
uncertainties in the reddening are too large for this approach to be any more certain than using the iron lines. 

Column 5 of Table 2 lists the derived MK classifications, which were found using line-depth ratios (e.g., 
Strassmeier and Fekel^^ or Gray and Johanson^^) and are not from the derived atmospheric parameters; however, 
it should be noted that in general the agreement between the physical parameters and the classification is quite 
good. 



Table 2. Derived stellar parameters 



Lensing Event 


Teff (K) 


log g (cm s ^) 


^ (km s-i) 


MK class 


Vrad (km s ^) 


97-BLG-45/47 


6000 


4.30 


2.48 


G5 V 


+59.98 ± 2.24 


97-BLG-56 


4480 


1.40 


1.84 


G8 III 


+177.0 ± 1.9 


98-BLG-6 


5200 


3.50 


2.28 


G8 IV 


-64.24 ± 1.68 


99-BLG-l 


5600 


4.65 


3.30 


G2 V 


+64.10 ± 1.78 


99-BLG-8 


4150 


0.95 


1.91 


K2 III 


+194.9 ± 1.7 


99-BLG-22 


5885 


4.30 


1.52 


FO V 


+37.35 ± 2.06 



4. RESULTS 
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In the Tables 3-6 we present an early look into the abundances of our program stars, with the caveat that 
some numbers are subject to change; i.e., elements that suffer from heavy hyperfine splitting (hfs) and/or 



isotopic shifts will likely be revised in ftiture work.^^ Currently, we correct for such effects when possible by 
correcting one or two lines with the atomic data of Prochaska^'' and applying an estimated correction to the 
total abundance. Abundances listed with a subscript (e.g., [Sc/Fe]i and [Sc/Fejn) are derived from both the 
neutral lines (I) and the singly ionized lines (II), and are listed separately. In all other cases the abundances 
come from lines of either ionization state. 

Tables 3-6 do not list uncertainties for abundances other than iron; typically errors from line-to-line scatter 
for a given element are ~ 0.2 dex but can vary quite a bit. This is especially true for silicon where the rms 
is more like 0.5 dex; however, [Si/Fe] is determined from lines with excitation potential around 6 e.V., while 
atmospheric parameters came from Fc I lines with x < 5 e.V.. On the other hand [Al/Fe] and [Na/Fe] often 
showed good agreement between the abundances derived from different lines (cr < 0.1 dex). 

We discuss the quality of the data for each star individually in Section 5. 



Table 3. Derived abundances 



Lensing Event 


[Fe/H] 


[Na/Fe] 


[Mg/Fe] 


[Al/Fe] 


[Si/Fe] 


[Ca/Fe] 


[Sc/Fe]i 


[Sc/Fe]„ 


97-BLG-45/47 


-0.28 ± 0.22 


+0.42 


+0.33 


+0.35 


+0.23 


+0.07 




-0.96 


97-BLG-56 


-1.13 ± 0.20 


-0.03 


+0.79 


+0.15 


+0.30 


+0.24 


+0.65 




98-BLG-6 


-0.22 ± 0.18 


+0.30 


+0.29 


+0.10 


+0.28 


-0.06 


+0.31 




99-BLG-l 


-0.01 ± 0.28 


+0.67 


+0.81 


+0.67 


+0.33 


+0.13 


+ 1.36 


-0.12 


99-BLG-8 


-1.12 ± 0.20 


+0.21 


+0.62 


+0.40 


+0.42 


+0.09 


+0.63 


+ 1.14 


99-BLG-22 


-0.35 ± 0.20 


+0.12 


+0.59 


+0.33 


+0.34 


+0.20 




-0.45 



Table 4. Derived abundances (continued) 



Lensing Event 


[Ti/Fc]i 


[Ti/Fc]ii 


[V/Fc]i 


[V/Fe]ii 


[Cr/Fe]i 


[Cr/Fe]ii 


[Mn/Fe] 


[Co/Fe] 


97-BLG-45/47 


+0.92 


+0.50 


+0.52 


+1.23 


+0.45 




-0.25 


+0.61 


97-BLG-56 


+0.67 


+0.83 


+0.61 




+0.45 


+0.14 


-0.25 


-0.01 


98-BLG-6 


+0.50 


+0.43 


+0.54 




+0.08 




+0.08 


+0.00 


99-BLG-l 


+1.19 


+0.39 


+1.33 


+1.11 


+0.55 




+0.81 


+0.46 


99-BLG-8 


+0.68 


+0.34 


+0.17 


+0.48 


+0.35 


+0.75 


-0.28 


-0.15 


99-BLG-22 


+0.70 


+0.47 


+0.52 


+1.47 


+0.01 




-0.27 


+0.60 



5. SUMMARY OF RESULTS 
5.1. 97-BLG-56 and 99-BLG-8 

We begin our discussion of the program stars with the cool giants. Both spectra have excellent signal-to- 
noise ratios, which comes as no surprise since these stars are already quite bright and being aided by lensing 
makes them that much brighter. They both contain the lowest [Fe/H] in the sample, raising the possibility 
of systematic errors in their temperatures affecting the iron abundance determinations. We doubt this is the 
case since our reduction techniques have been compared favorably with our so-called "standard stars" and our 
model parameters match well with the MK classifications derived via line-depth ratios. The radial velocities 
of these two stars are the highest in the sample, giving rise to speculation that they could be members of the 
Sagittarius dwarf spheroidal; however, using the spectral type to estimate the absolute magnitudes and some 
assumptions about reddening,"^ we derive distances to these stars of roughly 3.3 - 3.9 kpc. Compared with the 



Table 5. Derived abundances (continued) 



Lensing Event 


[Ni/Fe] 


[Zn/Fe] 


[Y/Fe]i 


[Y/Fe]ii 


[Zr/Fe] 


[Nb/Fe] 


[Mo/Fe] 


[Ru/Fe] 


[Ba/Fe] 


97-BLG-45/47 


+0.28 


+0.29 




+0.35 










+0.39 


97-BLG-56 


+0.25 




+0.84 


+0.42 


+1.14 




+1.14 


+0.92 


+0.57 


98-BLG-6 


+0.19 


+0.48 


+0.92 


+0.91 


+0.83 








+0.66 


99-BLG-l 


+0.31 


-0.22 


+1.43 


+0.96 


+1.76 




+1.31 


+ 1.94 


+0.01 


99-BLG-8 


+0.12 


-0.03 


+0.76 


+0.61 


+1.15 


+ 1.54 


+0.92 


+ 1.11 


+0.68 


99-BLG-22 


+0.23 


+0.03 




+0.05 










+0.26 



Table 6. Derived abundances (continued) 



Lensing Event 


[La/Fe] 


[Ce/Fe] 


[Pr/Fe] 


[Nd/Fe] 


[Eu/Fe] 


97-BLG-45/47 


+0.97 






+1.72 




97-BLG-56 


+0.81 


+1.27 


+0.90 


+0.99 


+0.53 


98-BLG-6 


+0.99 


+1.10 


+0.83 


+0.92 


+1.09 


99-BLG-l 


+1.47 


+1.55 




+1.31 


+1.46 


99-BLG-8 


+0.91 


+1.11 


+0.94 


+0.85 


+1.34 


99-BLG-22 




+2.39 


+0.89 


+0.89 





Sagittarius dwarf spheroidal, which has a distance of ^ 28 kpc,^'"^ these two stars clearly belong to a separate 
stellar population, even when accounting for errors in reddening or spectral type. 

The abundance trends in both stars are consistent: the low-Z (Z < 40) elements are enhanced by about 0.5 
dex while the high-Z elements are up by around 1 dex, indicating enhancements of s-process from possible AGB 

stars. The elements with estimated hfs corrections arc V, Mn, and Co for both stars, leaving the abundances 
from the high-Z elements perhaps overestimated; although, not all of these require hfs corrections. 

5.2. 98-BLG-6 

98-BLG-6 is roughly solar metallicity, particularly when the error bars are taken into account. The luminosity 
class is estimated from log g as either a main sequence turnoff star or a subgiant. The gravity does seem to 
be well determined given the excellent agreement between the abundances from the neutral and ionized atoms 
for the various species when available. The effective temperature/log g combination along with the spectral 
type from line ratios is consistent with this star being a subgiant (see, e.g., Thevenin and Jasniewicz^^), making 
98-BLG-6 the only one in this study. The radial velocity for this star differs from the previous two by 250 km s~^ 
and is the only one of the six that is negative. Estimated hfs corrections were done for V, Mn, and Co with 
corrections amounting to —0.65 dex for [Mn/Fe] and —0.45 dex for [Co/Fe]. Given these large corrections for 
the iron-peak elements, the high-Z elemental abundances should be viewed with caution. Finally, it should 
be noted that we are able to determine the abundances of 22 elements in a, V = 19.5 star only because of 
microlensing. Such detailed calculations are the exception at this magnitude, but should become trivial with a 
new age of 30-m and larger class telescopes, provided they are outfitted with the proper instrumentation and 
can function in the visible bands. 

5.3. 99-BLG-l 

According to our spectral typing, 99-BLG-l is a solar analog in the Galactic bulge. Our derived physical 
parameters are also consistent with solar-like parameters within their errors. Some lines from Sc II, V I, Mn I, 



and Co I have been corrected for hfs effects. Both [Sc/Fejn and [V/Fe]ii change very little as a result, which 
is surprising given the large overabundance of [V/Fc]n; however, there are 18 other V II lines that have yet to 
be corrected and it remains to be seen whether the one line we corrected at A5671 is anomalous or is actually 
indicative of large V enhancements. The high-Z elements as a whole show the largest enhancements of all the 
stars. Ruthenium is particularly enhanced; however, the abundance comes from only one line at A5636 and this 
element does have seven stable isotopes with a rather flat distribution in relative abundance in nature among 
them. Barium is the one anomalous high-Z element in this star with a scaled solar ratio; this despite it's sister 
clement. La, being enhanced by a factor of 30 relative to the sun. No hfs corrections have been applied to either 
of these elements. On the other hand, europium is enhanced and the high Eu/Ba ratio would indicate heavy 
r-process influences. Further analysis of this star should yield an interesting chemical history. 

5.4. 99-BLG-22 

This star is close to being a solar analog in temperature and gravity; however, the derived spectral type appears 
to make it too hot. The signal-to-noise ratio is quite good for this star (especially considering it has an unlensed 
V magnitude of 19.6), giving us confidence in our spectral analysis. Relatively little information is available for 
the high-Z elements, but once again, [Ba/Fe] is near solar while the heavier elements remain high. [Sc/Fe]n, 
[V/Fe]i, and [Mn/Fe]i have estimated corrections for hfs effects. Despite the quality of the spectrum, abundances 
of Zn, Y, Ce, Pr, and Nd come from only one or two lines each, which is atypical compared with the brighter 
stars. 

5.5. 97-BLG-45/47 

We conclude our summary with a discussion of 97-BLG-45/47. Like 99-BLG-22, the derived stellar parameters 
are at odds with the spectral type, and given the poor signal-to-noise ratio in the spectrum it is difficult to decide 
on a preferred set of parameters. The choices we adopted in Table 2 are found through spectroscopic means 
and are consistent with the parameters listed in Minniti et al.,^ making the star a GO dwarf. 97-BLG-45/47 is 
slightly underabundant in iron and shows a range of scatter in the remaining elements. [Nd/Fe] is particularly 
high, but the results do come from five lines and show only 0.34 dex scatter. The results for [Sc/Fe]ii, [V/Fe]i, 
and [Mn/Fc]i do include estimated hfs corrections. 

We revisited the A6708 Li I feature in the spectrum of 97-BLG-45/47 as first studied by Minniti et al..^ 
They measured an equivalent width of 47 mA and estimated an abundance of loge(Li) = 2.25 from published 
tables. We measure the equivalent width to be 37.3 mA, which when run in a full synthesis using the line list 
from Ford et al.,^'' results in loge(Li) = 2.25; i.e., good agreement with the earlier work. A question, however, 
arises as to the validity of the detection, given the noise in the spectrum and the fact that the center of the 
measured line is 0.12 A to the red of the strongest feature in the Li doublet, while the center of the Fe I line at 
A6705 is well matched in the data, as can be seen in Fig. 2. 

Finally, we note that our radial velocity is 30 km s~^ lower than that listed in Minniti et al.^ for uncertain 
reasons. Given the excellent agreement we obtained between the radial velocities of our "standards" observed 
on the same nights as 97-BLG-45/47 and their published values, we are quite confident in our determination 
listed in Table 2. 

6. CONCLUSIONS 

Microlensing has allowed us to derive a reasonably large and detailed set of abundances for a group of stars 
that have heretofore been beyond our capabilities. While some work on the details remains to be done before 
conclusions can be firmly drawn, we have been able to study abundance trends in stars as diverse as solar 
analogs and metal-poor red giants. Our early results show a trend of overabundance in nearly all the elements 
regardless of metallicity or spectral class. The two giants in the study also have the lowest metallicities and the 
highest radial velocities, and they are easily ruled out as members of the Sagittarius dwarf spheroidal. We have 
also re-analyzed the Li abundance in the star 97-BLG-45 and have found that we cannot conclusively rule for 
or against the presence of Li. All the spectra show some hint of Li A6708 and we will examine these in kind to 
determine abundances or upper limits. 
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Figure 2. The region around the Li I doublet at A6708 showing the spectrum of 97-BLG-45/47 (x's and dotted line) 
and the synthetic spectrum with loge(Li) = 2.25 and loge(Fe) = 7.42 (solid line). 

Our results arc likely to change for some of the elements where hfs or isotopic corrections are still necessary. 
When possible we corrected for hfs effects by correcting one line and applying a shift to the results from all the 
lines. Sometimes the shift was large (> 0.5 dex) and sometimes it was nonexistent. 

Thanks to a helping hand from Nature, we were able to get a glimpse at the functioning of Keck as a 
much larger telescope than its current 10-m diameter. Future 30-m and larger class telescopes will make 
such observations commonplace and allow us to learn truly the diverse histories of the Galactic bulge stellar 
population. 
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maps of the MACHO fields, and Tom Barlow for modifying his MAKEE routine to accept our data format 
and answering our many questions. This research has made use of the SIMBAD database, operated at CDS, 
Strasbourg, France. 
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